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D , Abstract. Electronic nematicity, proposed to exist in a number of transition 

metal materials, can have different microscopic origins. In particular, the 
anisotropic resistivity and meta-magnetic jumps observed in Sr3R.U207 are 

^ ' consistent with an earlier proposal that the isotropic-ncmatic transition is 

generically first order and accompanied by meta-magnctism when tuned by a 

S magnetic field. However, additional striking experimental features such as a 

non-Fermi liquid resistivity and critical thermodynamic behaviour imply the 
presence of an unidentified quantum critical point (QCP). Here we show that 
■"rt ' orbital degrees of freedom play an essential role in revealing a nematic QCP, 

C^ ' even though it is overshadowed by a nearby meta-ncmatic transition at low 

Q _ temperature. We further present a finite temperature phase diagram including 

f\ , the entropy landscape and discuss our findings in light of the phenomena observed 

in Sr3Ru207. 
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1. Introduction 

A variety of transition metal materials such as the cuprates [T| , Ru-oxides [2] , and Fc- 
pnictides [3] has been proposed to harbour an electronic nematic phase [4] . Electronic 
nematic phases are broadly characterized by the presence of spontaneously broken 
rotational symmetry and viewed as the quantum counterpart of nematic classical liquid 
crystal phases. The theoretical proposal of nematic quantum liquid crystals became 
more concrete when experiments on ultra-pure bilayer ruthenate (Sr 3 Ru 2 07) samples 
subjected to a magnetic field along the c-axis revealed an unusual phase characterized 
by a pronounced residual resistivity in place of a putative meta-magnetic quantum 
critical point (QCP) [5]. Interestingly, Sr3Ru2 07 was initially viewed as a prototype 
for the study of quantum phase transitions, exhibiting a striking non-Fermi liquid 
resistivity thought to originate from the putative magnetic field tuned QCP [6]. The 
unusual phase found in ultra-pure samples is delimited by two consecutive first order 
meta-magnetic transitions at low temperature and, remarkably, exhibits a significant 
in-plane magnetoresistive anisotropy when the external field is slightly tilted towards 
one of the in-plane crystal axes [2] . These observations strongly imply the formation 
of an anisotropic metallic, i.e. electronic nematic, phase in the bilayer ruthenate 
compound. 

Based at first on the two consecutive metamagnetic transitions, an electronic 
nematic phase was proposed and generic features of nematic phase formation were 
theoretically explored early on [SJ El [5] . It was found that the transition between the 
isotropic and nematic phase is generally first order, and that nematic order typically 
develops near a van Hove singularity (vHS) to avoid a Lifshitz transition. Varying 
the chemical potential, the nematic phase is bounded at low and high values by two 
isotropic phases, while the concomitant first order transitions lead to jumps in the 
electron density. When a magnetic field is applied (and the chemical potential is held 
fixed at, say, some low value), Zeeman coupling acts as a spin-dependent chemical 
potential term. Tuning the Zeeman field increases the volume of, e.g., the spin- up 
species which then enters the nematic phase near the vHS, resulting in a jump in the 
spin-up density at the isotropic-nematic transition. In contrast, the spin-down Fermi 
volume continuously decreases without passing a van Hove point. Such behavior gives 
rise to meta-magnetism, i.e. a sudden jump in the difference of up- to down-spin 
density at the isotropic-nematic transition [9]. Using more realistic band structures, 
microscopic models were suggested and the above conclusion was found to be robust 

nanu. 

Despite this mechanism of nematic phase formation, the origin of the critical 
signatures remains mysterious and recent thermodynamic data revealed additional 
complexity that requires new insight into the nematic theory |12j . For over a decade 
in temperature T, the resistivity behaves nearly perfectly linear [2], while the specific 
heat coefficient (C/T) varies as log(T) over the same range [321 [13] . This behaviour 
resembles quantum criticality observed in a variety of other correlated materials, 
suggesting the existence of an unidentified QCP related to the nematic phase. How 
can one reconcile the first order nature of the isotropic-nematic transition and the 
critical behavior associated with a second order transition? 

In this paper, we present a way to resolve this puzzle by introducing two nematic 
order parameters relevant for multi-orbital systems, and investigate as to how their 
interplay leads to interesting physics. While the current work is motivated by the 
bilayer ruthenate, our study on the isotropic-nematic transition is generally applicable 
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Figure 1. (a) Finite temperature phase diagram and contour plot of the entropy 
difference between the nematic and isotropic phase [S'(An) — S(0)]/T. The meta- 
ncmatic transition (marked by the arrow) and the QCP (marked by x) are shown, 
(b) The corresponding absolute entropy landscapes with (S'(An)) and without 
(5(0)) nematic ordering (for details see main text). 



to a t2 g -orbital system. We show that in the presence of moderate spin-orbit (SO) 
coupling the multi-orbital nature of the system is a key property for turning a nematic 
first order transition into a QCP. We also show that the nematic QCP is accompanied 
by a nearby jump in nematicity, dubbed here meta-nematic transition, which obscures 
the QCP at low temperature. This finding is in contrast to the current wisdom that 
a QCP is hidden under the nematic dome. 

Figure [T] (a) shows the finite temperature phase diagram of such a multi-orbital 
system, including the contour plot of the entropy difference between the nematic 
and the isotropic phase. The nematic QCP marked by x is found on the right 
low temperature phase boundary. Note also that the entropy inside the nematic 
phase near the meta-nematic transition indicated by the arrow is higher than in the 
isotropic phase. While the separation between the continuous and the meta-nematic 
transition depends on the underlying band structure and the SO coupling strength, 
thermodynamics and transport at finite temperature will be governed by the QCP 
aside from singular density of states (DOS) effects, which may be relevant for the 
phenomena observed in Sr3Ru2 07 [21 El H2J HH [HI [EJ . A more detailed discussion 
will be presented below. 
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2. Nematic QCP 

Theoretical progress has recently been made in developing microscopic models for 
nematic phase formation in Sr3Ru20y. One approach suggested that the nematic 
phase is a spontaneous Fermi surface (FS) distortion arising from a band with strong 
d xy character and a vHS near the Fermi level |llj . while another proposed that it 
originates from orbital ordering, i.e., the density difference between d xz and d yz 
orbitals, ignoring the d xy orbital [17U18J . Both scenarios lead to broken x-y symmetry. 
To distinguish both types of ordering in the present paper, wc call the former nematic 
and the latter orbital order. Although the underlying microscopic driving mechanisms 
are distinct, both approaches have noted the importance of SO coupling, which is 
in line with the results of angle resolved photocmission spectroscopy (ARPES) |19j . 
However, they failed to locate a QCP. Thus, given that all t2 g orbitals are coupled 
via SO interaction, we explore here a broader notion of x-y symmetry breaking in a 
multi-orbital system. 

To understand the origin of nematicity and its consequences, we start from a tight- 
binding model, which reproduces the FS of Sr 3 Ru207 [20l[2TJ[22] an d was introduced 
in [IT]. The model is based on the Ru 4d t2 g orbitals and includes moderate SO 
coupling Hso = ^XJ2 i L i • S^. Wc also incorporate a staggered lattice potential 

H s t = gJ2ka cr( c kI c k+Qcr + n - c wrtn Q = ( 7r , i") to account for the effect of the 
staggered rotation of the RuO@ octahedra, which doubles the real space unit cell 
[23l [TTl [24] . In principle, the lattice potential can be momentum dependent (e.g. 
g(k) ex cosk x + cosk y as in [23] ) but here we assume a constant value g as the results 
do not depend on its momentum structure. The tight-binding band structure then 
has the form 

*> = E «L ( A S .w ) Ct " (1) 

k,<7=t,4- 

/ el* 4 D + 10X -aX \ 

A ka = \ e™-iaX el* iX , G = gl 3X 3, (2) 

y — a\ —iX e£ y / 

where c£ CT = (c^J,c^, c^_ a , c^ T Qa , c^ Q(T ,a^ Q _ a ) consists of electron operators 
creating an electron with spin a =j, \ in one of the t2 g -derived orbitals a = 

yz,xz, xy. The orbital dispersions are given by e\^ xz = —2tiCosk y / x — 2t2Cosk x / y + (x, 
e k !/ = — %tz{cosk x + cosk y ) — 4t4COsk x cosk y — 2t^(cos(2k x ) + cos(2k y )) + /i, while the 
hybridization between the quasi-lD orbitals is e^ 13 = — 4<6sinfc.,;Sin£;j,. In the following 
the underlying band structure parameters are t\ = 0.5, ti = 0.05, £3 = 0.5, £4 = 0.1, 
£5 = —0.03, te = 0.05, g = 0.1 and A = 0.14 (if not stated otherwise), while the 
chemical potential (i serves as a tuning parameter. All energies throughout the paper 
are expressed in units of the intra-orbital nearest neighbour hopping integral 2t\ = 1. 
Since the nematic phase reported in Si'3Ru207 arises from a metallic state, 
it is reasonable to assume that long-range interactions are well screened, leaving 
moderately weak (compared with the bandwidth) on-site and nearest neighbour 
interactions. The microscopic interaction Hamiltonian is then 

H int = uj2 ?4< + \ E ™>-' + E va <<r n %> . (3) 

i,a i,a^/3 (i,j),ct 
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Figure 2. (I) & (II) represent nematic (A> 



,-! 



Ajj) and orbital (Aq) 



order parameters over the T-/x plane, respectively. Note the small magnitude of 
Ao- (HI) shows the qualitative changes of the FS in half of the original Brillouin 
zone at selected points in the phase diagram for T rs 0. In order of increasing 
/i these are taken just before and after the nematic transition at fi c i = —0.56 
(a & b, respectively), at fi = —0.535 (c), just before and after the meta-ncmatic 
jump at /i c 2 = —0.497 (d & e, respectively), and at the continuous transition at 
Mc3 = -0.47 (f). 



with nf a — c°^cf a the density operator for electrons in orbital a at site i and spin a. 
Here, U, U, and V a represent repulsive intra-orbital on-site, inter-orbital on-site, and 
intra-orbital nearest neighbour interactions, respectively. We assume that V xz (y z ) is 
finite along nearest neighbour x(y)-bonds, while V xy is finite along all four nearest 
neighbour bonds, and that V a = V. 

Different instabilities compete within Ho + Hi nt . Most natural candidates are spin 
density wave and ferromagnetic instabilities. However, it was shown in |17j that orbital 
ordering between d yz and d xz orbitals dominates when the inter-orbital interaction is 
significant between the two quasi-lD orbitals in the presence of SO coupling. This 
important observation, however, did not take into account the 2D d xy orbital, which 
dominates the 7 FS sheets of Sr3Ru207 as observed by ARPES [20]. In particular, the 
72 sheet possesses a vHS near the Fermi level implying that it is most susceptible to an 
instability. Reference [TT] indeed suggests that nematic ordering in the xy dominated 
bands is the leading instability when the nearest neighbour interaction is taken into 
account. Since these theories imply that ferromagnetic and antiferromagnetic spin 
density wave instabilities are suppressed by SO coupling, we focus here on the interplay 
of nematic and orbital ordering. 

Introducing nematic and orbital order parameters, A^ = A r_1 ^ k (cosfc x — 



cos k y )(n^ y a ) and A = N 1 J2 
Hamiltonian 



k.cr* 



-.XZ 
L kcr 



"k<T 



) , respectively, one arrives at the following 



H MF =H -V o J2[ A o( n 



XZ 
kcr 



"ka> 



i(Ao> 2 ; 
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Figure 3. Each panel represents the nematic order parameter for various SO 
coupling strengths A (dotted lines visualize discontinuities). Note that the first 
order transition at the right phase boundary becomes continuous for 0.07 < A < 
0.2. The results in figure [T] and [2] are based on A = 0.14. Very recent unpublished 
transport data |25| in fact bear a strong qualitative resemblance to this A = 0.14 
plot. The resistivity has a pronounced shoulder for fields higher than those of the 
two first-order transitions, and within this shoulder becomes anisotropic under 
the application of small in-plane magnetic fields. 



-V N ^[A N (cosfc x - 



k.rx 



cosfcj,)n^- 



2 (as) 2 ; 



(4) 



where the effective interactions are given by Vn = V and Vb = U/2 — 11/4 — V. 
Note that the on-site intra-orbital interaction U hinders orbital ordering (favouring 
magnetic ordering instead), but does not interfere with nematic ordering. On the 
other hand, the inter-orbital interaction U favours orbital ordering, while the nearest 
neighbour intra-orbital interaction V favours nematic over orbital order. 



The self-consistent solutions for nematic (Aj. 



k * 



A N ) and orbital order over 



the T — fi plane are shown in panels (I) and (II) of figure [2] Panel (III) shows the FS at 
selected points for zero temperature. As pointed out in figure [TJ nematic and orbital 
order develop continuously from the isotropic phase at ^ C 3 = —0.47, in contrast to 
the generic first order transitions encountered when either nematic pj [5J QT] or orbital 
[17] ordering is considered. 

For the present discussion we have set Vn = 0.8, while Vb(^C Vn) is chosen to 
be small or zero such that orbital order can only be induced by coupling to nematic 
order via SO interaction. However, orbital ordering can also be induced for larger Vq 
and smaller Vn values than the above choice (for example for Vb = 0.2 and Vn =0.6) 
but does not alter the main conclusion of the existence of the QCP and the nearby 
mcta-ncmatic jump. The reason for preferring the nematic to be the leading instability 
(Vn > Vb) is based on experimental constraints. We found that an independent orbital 
instability occurs for Vb > 0.55 (when Vn = 0) but that it is also accompanied by a 
drastic change in the overall FS topology, exhibiting various strongly ID-like open FS 
sheets inconsistent with the de Haas van Alphcn experimental findings for Sr3Ru207. 
In contrast, the nematic instability gives rise to smooth changes of the FS before and 
after the nematic window as shown in figure [2] (III), in agreement with the de Haas 
van Alphen results [211 E2] • Hence we consider the nematic instability as the relevant 
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instability in this study, yj 

To understand the effect of SO coupling, we investigate how the phase transitions 
are modified by changes in the SO coupling strength A. Figure [3] shows the nematic 
order parameter for different A (orbital order is not shown, but is small and disappears 
for A = 0). The nematic order parameter exhibits typical first order transitions for 
small A. For intermediate A, the first order transition at the higher critical chemical 
potential turns into a continuous transition preceded by a sudden meta-nematic jump. 
The change in nature from first order to second order is driven by the coupling between 
the two order parameters enabled by SO interaction and can be captured qualitatively 
by the following Ginzburg-Landau (GL) free energy analysis. 

3. GL free energy analysis 

A generalized GL free energy takes the form 

Tgl = a„A N + /? n A N + 7n A N + a A 2 + /3 A^ + 7o A^ 

+ <5A N A [1 + ??(An + A&)] + kA n A^ + . . . , (5) 

where we assume that the higher order terms denoted by the ellipsis are negligible. 
The GL free energy analysis is a phenomenological approach, and it is important to 
understand how each term is allowed by symmetry, in particular the S term. The 
nematic order parameter breaks only the x-y symmetry of the underlying lattice (or 
the rotational symmetry in the case of the Jellium model). Thus the GL free energy 
should be invariant under a sign change of the order parameter (e.g. An — > —An, 
corresponding to a tt/2 rotation), which restricts the GL expansion to only even powers 
in the nematic order parameter. Using the Jellium model of the electron gas, the 
GL free energy with nematic order was analyzed in |26) . where the isotropic-nematic 
transition is assumed to be continuous with a n changing sign at the transition, while 
/3 n stays positive. However, it was first pointed out in [8] using a single nematic order 
parameter that the free energy has a log-singularity due to a vHS (Lifshitz transition) 
leading to a first order isotropic-nematic transition in a simple square lattice model. 
An order parameter expansion shows that the /3 n coefficient becomes negative while 
a n and 7 n remain positive at the transition point [5]. Since then, turning the first 
order transition into second order has become an interesting problem. In the present 
study of multi-orbital systems such as the t2 g -orbital Sr3Ru207 compound, there are 
two x-y symmetry breaking order parameters as introduced in the previous section. 
Both orbital ordering (i.e. the density imbalance between d xz and d yz orbitals) and 
nematic ordering (defined on the d^y-orbital manifold) break the same x-y symmetry. 
In principle, the GL free energy contains a finite coupling term such as ~ An ■ Ao, 
as long as the Hamiltonian contains mixing between d xy and quasi- ID (d xz and d yz ) 
orbitals. However, due to the nature of the orbitals and the inversion symmetry of the 
square lattice about the x-y plane, finite hopping integrals between d xy and the quasi- 
1D orbitals are absent. Therefore, the dominant mixing between d xy and quasi-lD 
orbitals is SO coupling. 

Under the symmetry considerations discussed above, the 5-term is allowed when 
SO interaction is present, i.e. 5 oc A + 0(A 3 ). Starting with positive coefficients only, 
no ordering is initially present. However, as was shown previously for a single band, 

X While bilayer coupling may alter the role of orbital and nematic orders, the existence of a nematic 
QCP is not sensitive to the choice of leading and induced orders. 
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the nematic instability is typically of first order due to j3 n < and a n ^ n > [8]. 
Therefore, without SO coupling (S = 0) and a ,/3 > finite nematic order develops 
when f3 2 > Aa n j n , while orbital ordering is absent (Aq = 0). On the other hand, 
when 5^0 nematic order starts to develop continuously for S 2 > 4a n a , inducing 
also orbital ordering (with Ao = -^-An up to second order GL terms in Ao)- A 
jump in the order parameters furthermore occurs when, e.g., the quartic order 77-term 
becomes sufficiently negative. 

4. DOS, entropy and finite temperature phase diagram 

The finite temperature phase diagram is displayed in figure Q] (a). The phase boundary 
between the isotropic ground state and the nematic phase is continuous except near 
/i c i = —0.56, where a first order transition appears, extending up to T w 0.015 in 
temperature. Furthermore, the meta-nematic transition appearing at low temperature 
near fi C 2 = —0.497 disappears at a slightly higher temperature. Remarkably, the slope 
of the phase boundary near the QCP at /i C 3 = —0.47 is positive, indicating that the 
nematic phase has a higher entropy than the isotropic phase. 
To check this, we calculated the entropy 






S = ^ ^ [/(K) ln(/(K)) + (1 - f(K)) ln(l - /(££))] , (6) 

n k 

where f(E) is the Fermi function, which is also shown in figure [1] Over most of the 
region the entropy associated with the bare band structure is larger than the one for 
the nematic phase and peaks (or, rather, logarithmically diverges) at low temperature 
near \x = —0.54 due to the vHS in the DOS associated with the saddle points in the 
72 band as shown in figure Q] (b) . This is expected since previous work [5J [TU] , as 
discussed in the introduction, showed that the formation of the nematic phase helps 
to avoid a large DOS near the Fermi level, and hence a large entropy, by splitting a 
logarithmic DOS peak into two peaks located further away from the Fermi energy. 
One of the peaks is indeed shifted to lower /i values in the nematic phase, which is 
illustrated in figure 0] (see e.g. the 72 peak(s) in the panels for fx — —0.56 and —0.558 
near the first order transition at // c i)- Note also that the DOS peak deriving from 71 
band near the Fermi level clearly reacts to the formation of nematic order by splitting 
as well. 

There is another peak in the entropy landscape at the meta-nematic transition 
/i C 2 inside the nematic phase where the nematic entropy exceeds the entropy of 
the isotropic phase. This singular behaviour is unexpected and occurs since finite 
nematicity persists beyond fj. C 2 due to the coupling of nematic and orbital order as 
discussed in GL free energy analysis. This unusual feature is also responsible for the 
concave curvature of the phase boundary in the finite temperature phase diagram 
in figure Q] (a) . In the DOS the meta-nematic transition corresponds to a sudden 
changeover from one side of the Fermi level to the other of the lower of the two split 
72 van Hove peaks (cf. the sequence from /x = —0.498 to —0.492 in figure HJ, while 
only at the QCP (/x C 3 = —0.47) do the 7 vHS merge into one (respectively for 71 and 
72)- 
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Figure 4. Evolution of the DOS of all bands combined (black curves) and of the 
71 (red) and the 72 (blue) band (contributing most to the DOS near the Fermi 
level which is indicated by the vertical dashed line) with the chemical potential 
11 at T = (see also figures [T] and [2j . Note the splitting of the 7 band van Hove 
peaks near the Fermi level inside the nematic phase. 



5. Discussion and summary 

Quantitatively, the present model is based off the bilayer ruthenate compound, where 
the formation of a nematic phase is driven by a moderate magnetic field applied along 
the c-axis [2] . Although one of the tuning parameters of the present phase diagram 
is the chemical potential, let us attempt to relate our results to the experimental, 
magnetic field tuned phase diagram [12j . Within the weak coupling approach, the 
leading effect of a magnetic field along c-axis is to shift the location (in momentum and 
energy) of the relevant van Hove points in the 7 bands due to the Zeeman term |llj . In 
the absence of SO coupling, the Zeeman term largely acts as a uniform (in momentum 
space), spin-dependent chemical potential as discussed in section [2] In the presence of 
SO coupling, which mixes up- and down-spins, the effective chemical potential shift 
becomes momentum dependent. The shift is not uniform but determined by the spin 
composition of each band at a given momentum k. As discussed, the relevant 72 and 
71 vHS arise from a small region near (±7r,0) and (0, ±7r) in momentum space with 
a large d xy admixture. The effective shift of the ji S band (where s = ± denotes the 
pscudospin and i = 1,2) therefore is determined by the spin nature of the contributing 
d xy orbital (i.e. either cS. or cj/f). Hence we do not expect qualitative changes to our 
main conclusion, i.e. the existence of a QCP and a meta-nematic transition, when the 
tuning parameter is an external magnetic field oriented along the c-axis. 

However, one may expect that for a tilted field the interplay of SO interaction, 
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Zeeman term, and additional orbital couplings may affect the underlying band 
structure in a complicated manner. The effect of an arbitrary magnetic field on 
nematic, orbital and magnetic instabilities therefore is complex and requires a self- 
consistent theory with a number of different order parameters including independent 
spin-up and -down nematic and orbital orderings and orbital dependent magnetization 
channels, so that the present result has limited applicability. 

Moreover, the phase boundary at low /z and low T bends towards the nematic 
phase, opposite to the experimental observations. However, additional bilaycr effects 
possibly change the underlying band structure such that the curvature becomes convex 
due to DOS effects and may even give rise to an additional nematic QCP. Further 
studies of the interplay of nematicity, orbital ordering, bilayer coupling, magnetization, 
and SO coupling in the presence of a magnetic field (accounting for both orbital and 
spin effects) are highly desirable and will be presented elsewhere [27] . 

Naturally, critical fluctuations accompanying the nematic QCP give rise to a 
quantum critical fan at finite temperature. Considering DOS and band structure 
effects only, one can generally expect that thermodynamic quantities diverge at most 
logarithmically. In contrast, the experimental specific heat and entropy curves for 
Si'3Ru207 follow a power law oc H C /(H — H c ) as a function of magnetic field H when 
approaching the nematic phase, indicating that critical fluctuations are important 
[T2I 125] and substantiating our finding. In addition to entropy, critical non-Fermi 
liquid response was observed early on in resistivity measurements when approaching 
the nematic regime from finite temperature [6], exhibiting a ~ T a dependence with 
a ra 1.0 — 1.5. Theoretically, nematic quantum critical fluctuations in continuum 
and lattice systems have been investigated as well and shown to lead to non-Fermi 
liquid behaviour [2J3 1301 1311 US] ■ However, a critical fluctuation theory in the presence 
of a singular DOS and in particular the magneto-thermal critical behaviour in close 
proximity to a meta-nematic transition at low temperature is still missing. It is also 
worthwhile to note that magnetic susceptibility measurements will be most sensitive 
to the divergences at the first order transition and the meta-nematic transition inside 
the nematic phase, which obscures the QCP. Such a study lies beyond the scope of 
the present paper and will be addressed in the future. While the details of the current 
results are based on the band dispersion of the bilayer ruthenate, our GL analysis with 
multiple order parameters, describing both continuous and disguised first order (jump 
in the order parameters) transitions, is in principle applicable to a range of quantum 
phase transitions in itinerant multi-orbital materials. 

In summary, we have shown that SO coupling in a multi-orbital t2 g system can 
have a considerable impact on the formation of a nematic phase. Although the nematic 
instability arises from the d xy band, it also induces a density imbalance between 
d yz and d xz orbitals in the presence of SO interaction, which in turn changes the 
nematic phase transition from first to second order, such that the FS topology remains 
qualitatively unchanged across the transition. Under what circumstances a continuous 
or a first order phase boundary is more favourable in the presence of SO coupling likely 
depends on the detailed band structure and energetics as discussed in our analysis of 
the GL energy. Furthermore, since the continuous transition does not preempt a 
vHS, which in the present case happens via a nearby meta-nematic transition, the 
entropy diverges inside the nematic phase, causing a concave curvature of the phase 
boundary. The meta-nematic transition also hinders the observation of the QCP at low 
temperature, making further experimental analysis of the thermodynamics near the 
nematic phase boundaries in Sr3Ru2 07 necessary to confirm the present theoretical 
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proposal. 
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